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Abstract — In this paper we consider stationary Forced KdV 
equation with positive Forcing term. The supercritical solitary 
wave solutions of the stationary Forced KdV equation are 
obtained. In order to obtain the solutions the domain of the 
problem has been divided in to three parts; the left, the middle 
and the right parts. The solution on the left and the right parts 
are obtained by an analytical method. The solution on the 
middle part is expressed in the terms of Weierstarss elliptic 
function. We have designed computer programs using 
Mathematica to produce the solutions. The complete solution 
was found by matching the solutions of all the three parts. We 
have found out that there are four different solutions according 
to the values of the phase shift. Only one solution is positive. 
Further research can be carried out for negative forcing terms. 

Index Terms — Stationary Forced KdV Equation; 
Supercritical solution;, positive Forcin 


At] x +2a/]r] x +/h] xxx = fjx) -co<x<co (3) 

Equation (3) was first derived by Shen (1993). The solutions 
of equation (3) are categorized according to the value of 
A (Upstream flow velocity) as follows (Shen, 1993). 

(i) Supercritical stationary waves. They occur only when A is 
positive and sufficiently large. 

(ii) Subcritical stationary waves. They occur only when A is 
negative and sufficiently small. 

(iii) Unsteady periodic soliton radiation. This solution 
appears when |/l| is small. Such a solution is called the 
transcritical solution. 

II. Problem Statement 


I. INTRODUCTION 

The forced KdV equation is a first -order approximation of a 
long nonlinear surface wave in a channel flow of an in viscid 
fluid of constant density over a bump or a dent. This equation 
is normally written in the form (Gong and Shen, 1994) 

Vt + folx + ^ ar l 1 h + Ptfm = f x { x ) -°°<*<°° t> 0 ( 1 ) 

where A > 0, Ct < 0 and J3 < 0 are constant and f(x) is a 
given function called the forcing term which is differentiable 
and has a compact support i.e. it is nonzero only in a closed 
bounded set. Equation (1) was first derived by Akylas (1984) 
and in an asymptotically reduced result from Euler equations 
of fluid motion and the corresponding boundary conditions. 
The unknown Jj(x,t ) represents the first order elevation of 
the free surface of the fluid from it is equilibrium level. The 
forcing function f(x) is due to the bottom topography of the 
fluid domain (such as a bump or dent) or due to an external 
pressure on the free surface such as the wind stress on the 
surface of an ocean. In the absence of the forcing function 
term i.e. f(x) = 0 equation (1) becomes the familiar Korteweg 
deVries(KdV) equation 

Tj t + 2tj x + 2aijij x + [3r] )u = 0 -co<*<co t > 0 ( 2 ) 

when TJ t = 0 equation (1) becomes the equation of our 
concern, the stationary forced (sfKdV) equation. 
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In this paper we study supercritical solitary wave solutions of 
equation (3) with positive forcing terms in the rectangular 
bump. In this work we shall take the value of A equal to 3 
and the length of rectangular bump equals 2 in order to 
compare with the previous results. 

Solitary waves refer to any surface wave profile that dies out 
at infinity which means that the free surface elevation I] has 
the property 

77 — >0, rj x ^>0and 77 — > 0 as *—>±00 (4) 

Integrating equation (3) once with respect to the independent 
variable from -00 to x, we have 

Atj Aarj 1 + = f(x), -co<x<co, /t>0(5) 


with 77 — > 0, and ij x — >0 as x— >±co 

In this paper we shall study equation (5) and (6) with the 
following conditions 


3 1 

(i )a = — <0, B = — and A = 3. 

4 6 


(ii) /(*) = « 


1 , 


a ^ „ a 

< x < — 

2 2 


0 otherwise 


where a is positive constant representing the length of a 
rectangular bump. 
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In order to obtain the solution of equation (5) and (6) the 
domain of the problem has been divided in three parts. The 
left (x < -a/2), the middle (-a/2 < x < a/2) and the right parts (x 
> a/2). The complete solution was found by matching the 
solutions of all parts. Section 3 and Section 4 describes the 
method of finding solution on the left and the right part 
respectively. Section 5 describes the solution in the middle 
part. In Section 6 we determine the phase shift L g which 

depends on the forcing terms and the continuity conditions 
between each adjacent part. In Section 7 we combine the 
solutions obtained in Section 4 and Section 5. The 
conclusions of this paper are summarized in Section 8. 


III. Left-Side Solution 

In this Section we need to solve equations (5) and (6) in the 
region of zero forcing term by using analytical method. 
Meanwhile, the solution in the region of non zero forcing term 
will be discussed later in Section 5. 

Outside the region of the forcing term, equation (5) and (6) 
becomes 


Ar/ +ar/ 2 + pr/^ =0 ^ <x< A > 0 (7) 

77 — > 0 , and rj x — > 0 — > 0 as x—> ±00 ( 8 ) 


We need to find the solution in the negative region (x < -a/2) outside the rectangular bump. In this region 
equations (7) and (8) are written as 

Arj + ar/ 2 + J3rj" =0 — oo<x<— A>0 (9) 

77 — >0, and r/ x — >0— >0 as *—>-00 (10) 

where equation (10) is the half negative solitary wave conditions. We start by rearranging equation (9) in 
the form 


j3rj" = -Ar/ - ar/ 1 -co<x < — , A> 0 (11) 


/ drj 

where 77 = 

dx 


The procedure we adopt here in order to integrate equation (1 1) is that we write 1] as 


// / drj 

i] = 77 


dr/ 


( 12 ) 


Substituting equation (12) into equation (1 1) we obtain 


Pr/‘ dr/ 1 = (- Ar/ - arp^/dt/ (13) 


Equation (13) is simply a first order separable equation which can be integrated easily to get 

P ( / ) 2 A 2 & 3 , /-t 

-{d ) =-~V --77 +C 


(14) 


Multiplying both sides of equation! 14) by — we then obtain 

a 


3 P / /\ J/t 2 3 d 

TTVJ )=-— r J +~ c ( 15 ) 


3A 


2 a 


2 a 


p 

where — is positive constant and C is the integration constant. 

a 

By the half negative solitary wave conditions (10), equation (15) has a real double root, r o and r o is 
smaller than the third real root r. Then equation (15) can be written as 

M, 

2 a 


1 / x 2 

W) = {r-r/lr/-r 0 ),r 0 <r (16) 
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If we now let 


v = rj-r o 



then equation (16) can be written as 


r dv^ 

J 


= v 2 (si 



where s x = r - r o > 0 . 


( 17 ) 


Equation (17) can be integrated to yield 


Hence, the actual solution on the left side is 


V = s 1 sec /? 2 



r/(x) = 



(18) 


(19) 


where L g is the phase shift to be determined in Section 6. 
All graphs of these solutions are shown in Figure 1. In 
Section 6 we shall show how to calculate L a . For the purpose 

of illustration we now use the values of L n obtained in 

Section 6 to produce the appropriate solution of the sfKdV on 
the left 


The graphs shown in Figure 1 are the solution of the sfKdV 
equation on the left side when A = 3 and a = 2. The values of 
L 0 used for this purpose are -0.20752, -0.797065, -1.792481 
and -2.248243 for graph (a), (b), (c) and (d) respectively. 


side. 



Figure 1: Feft-Side Solution 


IV. The Right-Side Solution 

a 

In this section we determine the solution on the right side ( X > — ) that out side the rectangular bump. In 
this region equations (7) and (8) become 
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Xq + aq 2 + fiq" = 0 , ^<x<oo, 2 > 0 ( 20 ) 

77 — > 0 , ^ » 0 as x — » +oo (21) 

where equation (21) in normally termed as the positive half solitary wave conditions. 

The different between the left and the right region is simply on the solitary wave conditions. The same 
procedure as used for the left side can be applied to find the solution in the right and we thus obtain 

*7 W = - sec h 2 - L 1 1 (22) 

where the phase shift L x = —L 0 . All graphs of this solution are shown in Figure 2. 

The graphs shown in Figure 2 are the solution of the sfKdV equation on the right side when 2 = 3 and a = 
2. The values of L x used for this purpose are 0.20752, 0.797065, 1.792481 and 2.248243 for graph (a), 
(b), (c) and (d) respectively. 



find the solutions in the middle. This means that we need to 
V. The Middle Solution solve equations (5) and (6) in the region of the forcing terms. 

We will show that the solution in this region can expressed in 
At this stage, we have obtained the solutions at both sides terms of Weierstrass elliptic functions. The solutions we 
outside the region of the forcing terms discussed in Section 3 obtained ln Sectlon 3 ’ 4 - and 5 wl11 be combined to give the 
and Section 4, respectively. What is left now is to ful1 solutions. This will be presented in Section 7. 


M .a 

In this region, \X < — , the solution must satisfy 


Xq + aq 2 + fiq^ = 1 , X>0 ( 23 ) 
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We shall be dealing with equation (23). In order to solve this equation we need to deduce some conditions. 

j / a j a 

These conditions come from continuity of // and Tj at X = — — and .V = — . 


a 


Now at X = , we have 

2 


n 


' o' 


V 


f a} 
V 2y 


= rj a and rj =rj x (24) 


where 


3/1 7 2 I — ^ ( \ 

~2a S,Xh iJp {x+L - X <25) 




-/if a x 

y 


(26) 


At X = — we have 

2 


n 


6 

V 2y 
r 








n 




v z y 


r „\ 




or rj 


f 


v z y 


= -rj 




v z y 


(26*) 


Integrating equation (23) from -a/2 to x (< a/2), we find 


Equation (27) can be written as 


where 


By making the transformation 


Equation (28) is converted into 


^( 7 7 / ) 2 =T 1~~ T I 2 ~~ T J 3 +T ?o ( 27 ) 


{d ' )' = b x j f + b 2 ?f + btf + Z? 4 (28) 


b, = b = -— b, = — and b d = b-Ji n 

1 3/3 2 /3 3 /3 4 3 ° 


rj = c,t/ + c 2 


= ^(qM + c 2 ) 3 +b 2 {c l u + c 2 ) 2 +b 2 {c l u +c 2 ) +b A (29) 

Dividing both sides of equation (29) by (cj )~ we obtain 
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f { i \ 

b 2 


(u) =(Z? 1 c 1 )m 3 + 

(W +b 2 c 2 2 + b 3 c 2 +b 4 )— 


3b, 

V \ 3 V 


+ b n 


\ ( ( , 2 A 

3b b 


2 , 
U + 


J 


9b ' 

V v yu \ 


+ 2 b. 


( i A 

b 2 


3b j 


+ b % 


v ; 


u 

— + 
Cl 


(30) 


Equation (30) is now written in the a simplified form as 


u' = Au i -g 2 u-g 3 (31) 


where 


c i , ’ C 2 ’ 8 


1 \ 


3b, 


and (W + V,^ 4 ) 


In equation (31), g 2 is constant and g 3 is a function of L o for given A, CC, fi and a. The general 
solution of equation (31) can be expressed in term of Weierstrass elliptic function (Gong, 1994) 

u(x) = fc?(x + T,g 2 ,g 3 ) 

where T is a constant. 


Thus, in the region of the non-zero forcing term ( |x| < — ) the supercritical solitary wave solution can be 
written as 


Tj(x) = c l p(x + T,g 2 , g 3 ) + c 2 


(32) 


Equations (25) and (26) are now respectively reduced to 


=c 1 $K-- + T, g 2 ,g 3 ) + c 2 (33) 
and 


t h=c l ^'(-^ + T,g 2 ,g 3 ) (34) 


From Gong and Shen (1993), we find the identity 

, ,rr v 1 \fp l (x,g 2 ,g 3 )-fp'(T,g2,g i )\ 2 r , , , trr 

(x + T,g 2 ,g 3 ) = -\ — — -&Kx,g 2 ,g 3 )-p(T,g 2 ,g 3 ) 

4 [ P(x,g 2 ,g 3 )-p{T,g 2 ,g 3 ) ) 

Substituting equations (33) and (34) into equation (35), we then obtain 


(x + T,g 2 ,g 3 ) = ^ 


^'(x,g 2 ,g 3 )~ 


Vx 


P(x,g 2 ,g 3 )-( T1 ° Cz ) 


- -^(x,g 2 ,g 3 )~— — — (36) 


(35) 
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We now observe that in equation (36), the right hand side is independent of constant T. We then have 


(x + T,g 2 ,g 3 ) 


1 j c 1 p'(x,g 2 ,g 3 )-r/ ] 
4{p(x,g 2 ,g 2 )-Tj 0 +c 2 


2 


~P(x,g 2 ,g 3 ) 


n o -c 2 

c i 


(37) 


Using Mathematica program Substituting equation(37) into equation (32), we get the solution on the 
middle side. All graphs of these solutions are shown in Figure 3. 

Figure 3 shows the solutions of the sfKdV in the middle, when A = 3 and a = 2. The values of L a are 
-0.20752, -0.792481, -1.792481 and -2.248243 for graph (a), (b), (c) and (d) respectively. 



VI. The Phase Shift L o 


In this section we shall determine the phase shift L 0 which depends on the forcing terms and the continuity 
conditions between each adjacent part. 

Substituting x = a/2 and x = -a/2 into equation (32) we get 


V 






(a A 

PP ~ + T ^2^3 


+ 


r a 7 


and V ~- + T,g 2 + g 3 
V A 


+ c 2 . 


Substituting the above expression in equation (26*) we obtain 

P 

From equations (33), (34), (37) and (38) we can obtain 


r a A 

~ + T,g 2 ,g 3 


(a A 

P -~ + T,g 2 + g 3 
V 1 


= 0 (38) 
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bU,l.) 


c, j C]P'(a,g;,g 3 )-7i 
4 [cMa,g 2 ,g 3 )-7io+ c 2 


2 


-Cip(a,g 2 ,g 3 ) + 2 c 2 


2rj o =0 (39) 


Equation (39) determines the phase shift L a . 

Equation (39) have a solitary wave solution only when B(A,L (l )=0 (Gong, 1994). 
Hence, B(A, L a ) = 0 is the condition used to calculate the phase shift L a . By using Mathematica 
program the function B(A, L a ) was plotted versus trial L 0 . The intersection point with L n -axis was 
found as the values of L a .(see Figure 4) 

This curve has four intersection point with L a -axis, i.e. B(A, L a ) has four zero, L 0 = —0.20752, 
-0.797065,-1.792481 and -2.248243 that means we have four supercritical solitary wave 
solutions. 



Figure 4: The Curve of B( A, L ( . ) 


versus 


VII. Complete Solutions 


solutions of the sfKdV equation. This means that we need to 
combine the solutions in both side with the middle ones. 


We have obtained the solutions of the sfKdV equation at both 
sides and also in the middle as discussed in Section 3, 4, and 
5, respectively. In this Section we shall give the complete 


We rewrite the solutions outlined in Section 3 as 


i A x )= (40) 


for the left side, and 


for the right side. 


v(x)= ~^ sech2 J^( x - L ')’ (41) 
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T?(x) = c l p(x + T,g 2 , g 3 ) + c 2 (42). 


By matching equations (40), (41) and (42) we get 


7700 = 


32 2 2 a 

seen ( x-L ) , -oo<x< — ; 

2a \ 4/3 ° 2 

. „ a a 

c 1 {p(x + T,g 2 ,g 3 ) + c 2 - — <■*<—; 


(43) 


32 2 2 a 

seen (x-L,), — <x<co. 

2 a \ 4/T 1 2 


Equation (43) is the complete solution of sfKdV equation. 

The graphs shown in Figure 5 are the complete solutions of the sKdV equation when 2 = 3 and a = 2. The 
values of L o are -0.20752, -0.797065, -1.792481 and -2.248243 for graph (a), (b), (c) and (d) 
respectively 



VIII. Conclusion 


Figure 5: The Complete Solution 

four different solutions according to the value of the phase 
shift L o . Only one solution is positive the three others 


In this paper we have studied the solutions of the stationary 
Korteweg-deVries (sfKdV) equation. In order to obtain the 
solution of the sfKdV equation, the domain of the problem 
has been divided into three parts, the left, the middle and the 
right parts. The solutions on the left and right parts were 


solutions have positive and negative values.. Further research 
can be carried out for negative forcing terms 
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